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ABSTRACT 

We examine the properties of a previously discovered population of globular clusters in the heart of the rich galaxy cluster Abell 
1185 that might be intergalactic in nature. Deep images obtained with the Advanced Camera for Suri'eys (ACS) aboard Hubble 
Space Telescope (HST) confirm the presence of ~ 1300 globular clusters brighter than Ifshw - 27.3 mag in a field devoid of any 
large galaxies. The luminosities and colors of these objects are found to be similar to those of metal-poor globular clusters observed in 
many galaxies to date. Although a significant fraction of the detected globular clusters undoubtedly reside in the outer halos of galaxies 
adjacent to this field, detailed modeling of their distribution suggests that the majority of these objects are likely to be intergalactic, in 
the sense that they are not gravitationally bound to any individual galaxy. We conclude that the true nature and origin of the globular 
cluster population in the core of All 85 - galactic residents or intergalactic wanderers - remains uncertain, and suggest how future 
observation could resolve this ambiguity. 

Key words, galaxies: elliptical and lenticular, cD - galaxies: formation ~ galaxies: interactions - galaxies: clusters: individual: Abell 
1 185 - globular clusters: general 



1. Introduction 

Globular clusters are ubiquitous in galaxies, found in all but the 
smallest dwarfs (see, e.g.. Figure 3 of Peng et al. 120081) . For this 
reason they have long played a key role in astronomers' quest to 
understand the origin and evolution of galaxies (see West et al. 
12004, and Brodie & S trader 2Q06. for reviews). 

However, a growing body of observational evidence also 
hints at the existence of a population of intergalactic globular 
clusters (hereafter IGCs) that reside outside of galaxies. Half a 
century ago, van den Bergh ( 119581 ) speculated that "about one 
third of all the globular clusters are of the intergalactic type" and 
numerous studies since then have lent support to the idea that 
at least some globular clusters are not gravitationally bound to 
individual galaxies (e.g., Muzzio et al. 1984, White 1987; West 
et al. [19951 Ba ssino e t. aL l2003l Hilker^Ml Jordan et aL l2003l 
WilUams et al.l2007l Schub erth et aL l2008l Coenda, Muriel & 
Donz elli l2009l Gregg et al. l2009l Lee et al. 12010: Peng et al. 

Eon). 



* Based on observations with the NASA/ESA Hubble Space 
Telescope obtained at the Space Telescope Science Institute (STScI), 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. 



Little is known at present about the nature or origin of IGCs. 
One possibility is that they are normal globular clusters that 
were stripped from galaxies during gravitational interactions, 
the detritus of cosmic evolution. Globular clusters, being among 
the densest stellar systems, are likely to survive the partial or 
even complete disruption of their parent galaxy, resulting in an 
ever-growing population of orphaned star clusters that roam the 
space between galaxies. Alternatively, some globular clusters 
might have been born in intergalactic space without ever resid- 
ing in galaxies, formed perhaps from the gravitational instability 
of primordial density fluctuations (e.g., Peebles & Dicke [T9681 
Peebles 1984; West 1993; Yahagi & Bekki 2005; Griff'en et al. 
2010 or thermal instabilities in cooling flows and other gas-rich 
environments (e.g., Fabian et al. 119841 Fall & Rees 119851 Cen 
i2001 ; Griff^en etaL|20l0l). 

Whatever their origin, IGCs are likely to be most abun- 
dant in the extreme environments of galaxy clusters, where tidal 
forces and collisions between galaxies can pluck globular clus- 
ters from their halos or the dense local conditions might have 
allowed the birth of massive star clusters from intergalactic gas. 
Unfortunately, searches for IGCs in galaxy clusters are hindered 
by the fact that most have a giant elliptical at their center where 
the density of IGCs is expected to be greatest, making it difficult 
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Fig. 1. A composite BVR image of the rich cluster Al 185 obtained with the CFH12k camera on CHFT (see Andreon et al. l2QQ7l t 
showing the locations of our ACS fields. The field of view is ~ 20' x 16', corresponding to ~ 0.8 Mpc x 0.64 Mpc at the distance 
of Al 185. North is up and East is to the left. Our primary ACS field is located just below the brightest cluster galaxy, NGC 3550. 
A second ACS field located ~ 9' (~ 350 kpc) to the northeast was used as a control field, as discussed in the text. Several of the 
brightest galaxies near the ACS field are identified, including the spectacular interacting galaxy pair Arp 105. X-ray contours from 
Mahdavi et al. ll996l based on Einstein observations are overlaid on the image. 



to distinguish genuine IGCs from the galaxy's intrinsic globular 
cluster population. 

An exception is Abell 1185 (z = 0.033), a rarity among 
low-redshift galaxy clusters because its brightest member galaxy 
(NGC 3550) is offset by ~ 150 kpc from the peak of the x-ray gas 
emission that presumably marks the centre of the cluster mass 
distribution (Mahdavi et al. 1996). Hence, Al 185 offers an ideal 
opportunity to search directly for IGCs in the core of a relatively 
nearby cluster of galaxies. Tentative evidence of such a popu- 
lation was presented by Jordan et al. (2003), who detected 99 
candidate IGCs in HST/WFPC2 images of a field in the center 
of All 85. However, their data suffered from incompleteness at 
magnitudes fainter than Ifm4w - 25.5 mag, allowing only the 
brightest ~ 10% of the expected globular cluster population to 
be detected at the distance of Al 185. 

In this paper, we analyze deep F555W and F814W ACS im- 
ages of the core of Al 185. We confirm the presence of a substan- 
tial population of globular clusters in this field and extend the re- 
sults of Jordan et al. (2003) by presenting the first observations 
of the luminosity function and color distribution of these ob- 



jects. The observations and data reductions are described in §2. 
Results from our analysis are presented in §3, and their implica- 
tions are discussed in §4. We assume a standard cosmology with 
Ho - 72 km s"' Mpc"' and a distance modulus m - M ^ 35.7 
mag to A1185. At this distance, an angle of 1' subtends a pro- 
jected linear scale of 65 kpc, and the 0.05" pixel scale of the 
ACS/WFC coiTesponds to a physical scale of 32 pc. 

2. Observations and Data Reductions 

Deep images of a field near the location of the peak x-ray emis- 
sion in A1185 (RA = llh 10m 46s, DEC = +28°43'57") were 
obtained using the ACS Wide Field Channel (WFC) onboard 
HST as part of program GO-9438, which was allocated 15 or- 
bits of time (see Figures [T] and |2]i. Two filters were used, F555W 
and F814W (roughly equivalent to Johnson V and Cousins /), 
with total exposure times of 15,540 s and 23,160 s, respectively. 
The longer integration time of the If^i4w images was chosen to 
ensure complete detection of the brightest half of the globular 
cluster luminosity function (GCLF) at the distance of A1185, 
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Fig. 2. ACS image of the central field in Abell 1 185, produced using the V^fsssiv and Ifsi4w images and averaging them to produce 
the third color. The brightest cluster galaxy, NGC 3550, is located just off the upper right corner of the image. The field of view is 
3.4' X 3.4', corresponding to a spatial scale of ~ 220 kpc on a side. North is up and east is to the left. The white box indicates the 
region shown in Fig.|5] 



which is assumed to have a Gaussian distribution with a peak 
at Mv - -7.4 mag (e.g., Harris 2001; Jordan et al. 2007) and 
hence should occur at an apparent magnitude Ipniiw - 27.3 mag 
for typical globular cluster colors (V-I ^1). The relatively shal- 
lower limit of the Vf^s^w images was intended to provide color 
information for brighter globular clusters only. With a scale of 
0.049" per pixel, the images have a 202" x 202" field of view 
corresponding to ~ 130 kpc on a side at the distance of A1185 
(for an assumed distance modulus m - M = 35.7). Individual 
exposures were not dithered because the small gap between the 
ACS chips does not adversely affect our ability to detect globular 
clusters over most of the field. 

Images were processed using the Apsis ACS data pipeline 
(Blakeslee 120031) . which produces final sky-subtracted images 
in each filter after distortion correction, image alignment, cos- 
mic ray removal and image combination using sub-pixel driz- 
zling (Fruchter & Hook 2002 ) with a final pixel scale of 0.05". A 
combined Vf555w + Ifshw image is shown in Fig.|2] Weight im- 
ages were also constructed for each output science image; these 



give the total uncertainty for each pixel based on instrumental 
and photon shot noise (see Jordan et al. 2004). These weight im- 
ages provide a map of the nns variation in the final image that 
was used for object detection as described below. 

2. 1 . Object detection and classification 

Object detection was performed using the Source Extractor 
(SExtractor) software package (Bertin & Arnouts |l996b to pro- 
duce a catalog of candidate globular clusters. At the distance 
of A1185, globular clusters are unresolved even with ACS and 
hence will appear as faint point sources (a typical globular clus- 
ter with half-light radius r/, ^ 3 pc would have an angular size 
much less than one ACS/WFC pixel, however this is diffused over 
several pixels because of the ACS PSF which has a FWHM ~ 2 
pixels). 

Objects were selected using a detection threshold of five con- 
tiguous pixels each at least 1.5cr above background after con- 
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Fig. 3. Top: SExtractor classes of all detected objects in the core 
of A1185. Bottom: SExtractor classes of artificial point sources 
that were added to the field as described in Sect. 2.3. Our adopted 
criterion of CLASS > 0.8 for globular cluster candidates does an 
excellent job of identifying point sources down to the faintest 
magnitudes. 



volving the image with a Gaussian filter gauss_2.0_5 x 5. The 
nns maps produced by Apsis were used to weight each pixel 
based on its noise level by setting the WEIGHT- TYPE parameter 
to MAP_RMS when running SExtractor. The local background was 
estimated using a mesh size of 40x40 square pixels. These detec- 
tion parameters were chosen after experimentation as providing 
the best compromise between sensitivity and reliability. Object 
detection was done independently for the F555W and F814W 
images and the resulting object catalogs were then matched us- 
ing a 0.1" matching radius. 

SExtractor assigns each detected object a classification pa- 
rameter, CLASS, that ranges from ~ for galaxies and other ex- 
tended objects to ~ 1 for point sources. CLASS values are based 
on a neural network algorithm trained on simulated images (for 
details see Bertin & Arnouts 1996 ). We consider point source 
candidates to have CLASS > 0.8. Fig.[3]shows that this criterion 
does an excellent job of discriminating globular clusters from 
galaxies and other non-point sources even at the faintest magni- 
tudes. Experiments in which artificial point sources of different 
magnitudes were added at random locations in the ACS field (see 
Sect. 2.3 for details) show that fully 98% of all input objects de- 
tected by SExtractor down to /fgMiv = 28 mag were assigned 
CLASS > 0.8. 



2.2. Photometric calibration 

Because globular clusters are unresolved at the distance of 
A1185, we used the ALLSTAR routine in the DAOPHOT soft- 
ware package (Stetson 1987 ) to perform photometry of detected 
candidates based on PSF fitting, which provides more reliable 
magnitudes than the default SExtractor magnitudes. 



To do this, the coordinates of SExtractor-detected objects 
were fed into DAOPHOT/ALLSTAR and instrumental magni- 
tudes were then measured at each location using PSF models de- 
termined empirically from stars on the final F555W and F814W 
images. We followed the procedures of Sirianni et al. (2005) for 
aperture corrections and converted instrumental magnitudes to 
the VEGAMAG system using the following zeropoints: 



555W = -2.5 log fsssw + 25.724 
/814W = -2.5 log /8I4W + 25.501 



(la) 
(lb) 



where fsssw and /suw are the integrated fluxes in units of elec- 
trons per second in each filter Reddening for the A1185 field 
was taken to be E(B - V) = 0.029 mag based on the DIRBE 
maps of Schlegel et al. (1998) and extinction ratios were taken 
from Sirianni et al. (i2005i . 

2.3. Completeness limits 

Object detection and classification become increasingly unreli- 
able at faint magnitudes. To assess our completeness limits, ar- 
tificial stars of known magnitudes were added at random loca- 
tions in the final ACS images and then the same SExtractor and 
DAOPHOT procedures described above were used to determine 
the fraction of stars recovered at each magnitude. These experi- 
ments were performed by adding 100 stars at random to the ACS 
field using the DAOPHOT task ADDSTAR. SExtractor was then 
run with the same detection parameters and the fraction of in- 
put stars that were recovered was recorded. This process was re- 
peated many times to generate the completeness function shown 
in Fig.|4]for the F814W observations. 

These experiments show that object detection and classi- 
fication are essentially 100% complete down to the expected 
turnover magnitude of the GCLF at Ifuhw - 27.3 mag, and 
declines rapidly at fainter magnitudes. To avoid uncertainties 
caused by incompleteness at fainter magnitudes, we adopted a 
completeness limit at If&i4w - 27.3 mag and consider only ob- 
jects brighter than this in subsequent analysis. 



3. Results 

Globular clusters detected in a small portion of the Al 185 core 
field are shown in Fig. |5] A total of 2285 such objects se- 
lected solely on the basis of their SExtractor classifications were 
found in the F814W image, of which 1422 are brighter than our 
adopted completeness limit //?8i4w - 27.3 mag. The projected 
spatial distribution of such objects is plotted in Fig.|6] Matched 
detections of 1 124 objects were found in the F555W image. 

It is clear from Fig. |6] that the two-dimensional distribu- 
tion of IGC candidates is not uniform. A prominent enhance- 
ment in the direction of NGC 3550 is easily seen, as are sev- 
eral other concentrations (e.g., around x ^ 500, y - 300 associ- 
ated with the galaxy MCG+05-27-004 that lies just outside the 
ACS field). Hence we are undoubtedly seeing some contamina- 
tion from globular clusters in the halos of galaxies, which we 
address below. 



3.1. The number of IGCs in A1185 

Jordan et al. (2003) detected 99 IGC candidates in Al 185 based 
on shallower WFPC2 observations. If intergalactic globular clus- 
ters have the same luminosity function as their counterparts in 
galaxies, then scaling from the WFPC2 results to the depth and 
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Fig. 4. Completeness of point source detection as a function of 
apparent magnitude for the F814W image of the central region 
in A 11 85. The red line shows the expected turnover magnitude 
of the globular cluster luminosity function at /f8I4w - 27.3 mag. 

area of our ACS field (assuming a constant density of objects 
across the field) we would expect to find approximately 2000 to 
4000 globular clusters in this region, with about half this num- 
ber detected down to our adopted 100% completeness limit at 
Ifsi4w - 27.3 mag. The detection of 1422 objects brighter than 
this magnitude limit is consistent with this estimate. 

However, it is inevitable that some IGC candidates are actu- 
ally foreground stars, unresolved background galaxies, and glob- 
ular clusters in the halos of galaxies near the ACS field mas- 
querading as IGCs. Confirming the presence of IGCs in A1185 
and estimating their number requires correcting for these various 
sources of contamination. 

To assess the contribution from stars and unresolved galax- 
ies along the line of sight, four fields observed with the same in- 
strument, filter and comparable integration times were processed 
and analyzed in an identical manner as the Al 185 core field. The 
selected fields, which are all at similar Galactic latitudes as the 
Al 185 field, are devoid of bright galaxies and can therefore pro- 
vide unbiased estimates of the number of contaminating point 
sources at difi'erent magnitudes. The four fields used were: 

- A second field in Al 185 (see Fig.[T]l was observed as part of 
our programme GO-10277. This field is located ~ 350 kpc 
from the primary field and hence can provide a useful con- 
straint on the spatial extent of the IGC distribution within 
A1185. ACS observations with integration times of 17,816 
seconds were obtained for both the F555W and F814W fil- 
ters. 

- F814W images of a field observed to study high-redshift Ly- 
a emitters (GO-941 1) were extracted from the HST archive. 
This field is at Galactic latitude b - 68°, almost identical 
to that of A1185. A subset of 16 images were combined to 
attain a total exposure time of 24,060 seconds. 

- Two additional fields from the Gemini Deep Deep Survey 
(GO-9760) were obtained from the HST archive. These are 
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Fig. 5. Intergalactic globular cluster candidates (in boxes) found 
in a small portion of the ACS field in the core of Abell 1 185. The 
field of view is ~ 0.3' x 0.5'. 

the SA02-1 and SA02-2 fields, which are located well away 
from the Galactic plane at b - -61°. The total integration 
time for both fields was 16,345 seconds with the F814W fil- 
ter 

Figure 7] shows the completeness functions for each of the 
control fie ds determined from artificial star experiments. All 
four show similar completeness limits, with nearly 100% of ob- 
jects detected down to Ifsi4w - 27.3 mag, which is not surpris- 
ing given their similar integration times. 

In Figure |8] we plot the number of detected point sources 
as a function of If^i^w magnitude for all fields. The agreement 
between the four control fields is remarkably good, which sug- 
gests that they provide a reliable estimate of the number of con- 
taminating foreground and background point sources along the 
line of sight to A1185. Furthermore, the fact that the counts in 
the second (outer) A1185 field are statistically indistinguishable 
from those in the three other control fields indicates that globular 
clusters are not present in significant numbers in this region, and 
hence if an IGC population is present in All 85 then it must be 
concentrated towards the innermost regions. 

Averaging results from the four control fields, we conclude 
that the number of contaminating foreground stars and back- 



6 



Michael J. West et al.: The Globular Cluster Systems of Abell 11 85 



4000 



3000 

I — I 
m 
X 

^ 2000 
1000 











. . ■ ■■■ V ■ 7 ■ 





















1000 



2000 
X [pixel] 



3000 



4000 



Fig. 6. The distribution of the 1422 intergalactic globular cluster 
candidates brighter than If&hw - 27.3 mag in the A1185 core 
field. The field of view is 3.4' x 3.4', with north up and east to 
the left. Enhancements in the number of objects can be seen in a 
few regions, especially in the direction of the bright galaxy NGC 
3550. For comparison, the location of the smaller and shallower 
WFPC2 field studied by Jordan et al. 2003 is also shown. 
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Fig. 8. The number of detected point sources in the Al 185 core 
field and the four control fields. Counts have been corrected for 
incompleteness at each magnitude. A large excess of objects is 
clearly seen at all magnitudes in the core field compared to the 
control fields, consistent with the presence of a substantial pop- 
ulation of globular clusters. 
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Fig. 7. Completeness of point source detection as a function of 
Ifs\4w magnitude for the four control fields. The red line shows 
the expected turnover magnitude of the globular cluster luminos- 
ity function at IfS]4w - 27.3 mag. Detection of point sources is 
nearly 100% complete down to this magnitude in all four control 
fields. The completeness fraction at faint magnitudes is greatest 
for the GO-9411 field, which has the longest total integration 
time. 



ground galaxies is ^ 142+12 down to a magnitude Ifsuw — 27.3 
mag. The number of objects detected in the A1185 core field 
is an order of magnitude greater than this, a total of ~ 1300 
globular cluster candidates. Likewise, Figure [9] shows the color- 
magnitude relation for point sources in the A1185 core field 
compared with that of similar objects in the outer field of Al 185. 
If we consider only those objects with colors in the range 0.7 < 
V -I < 1.5 expected for globular clusters, then there are 977 ob- 
jects in the central field compared to 32 in the outer field, strong 
evidence of the presence of a large population of globular clus- 
ters in the core of Al 185. Based on our detection of 1422 glob- 
ular cluster candidates brighter than If&hw - 27.3 mag, after 
subtracting the expected number of contaminating background 
objects and doubling the resulting number to include the faint 
half of the GCLF, we estimate that ~ 2560 globular clusters in 
total are present in the ACS field. 

To estimate the fraction of the detected globular clusters that 
are associated with galaxies, we performed detailed simulations 
of the expected distribution of globular clusters belonging to all 
galaxies brighter than V ^ 18 {My - -18) within 3' of the 
centre of the A1185 core field (fainter or more remote galax- 
ies are unlikely to contribute contaminating globular clusters in 
significant numbers). These 14 galaxies and their properties are 
listed in Table 1 along with their apparent V-band magnitudes 
measured from the CFH12k images of Andreon et al. 2007 us- 
ing SExtractor and their redshifts when known. Although two of 
these galaxies do not have measured redshifts, we have assumed 
that all are members of A1185. The following procedures and 
input parameters were used to generate the simulations: 

- It is well established that the number of globular clusters be- 
longing to a galaxy scales with its luminosity. A standard 
measure is the specific globular cluster frequency, S n, de- 
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Fig. 9. Color-magnitude diagram for globular-cluster-like ob- 
jects in A1185 core field. Black points are objects in the core 
field and red points are from the outer field indicated in Fig. 
[T] The prominent swath of objects with colors 0.7 < Vfsssw - 
Ifs\4w ^ 1-5 shows the detection of a substantial number of 
globular clusters in the core field. Within this color range there 
are 977 objects in the central field compared to only 32 in the 
outer field, a factor of more than 30 enhancement in the num- 
ber of objects. Mj denotes the absolute magnitude of detected 
objects assuming a distance modulus m - M - 35.7. 



fined as 



Sn = Ngr 10 



-0.4(Mv + 15) 



(2) 



where Ngc is the total number of globular clusters and My is 
the galaxy's absolute V-band magnitude (Harris & van den 
Bergh 119811) . Blakeslee, Tonry & Metzger 1997 measured 
the specific frequency of the dominant elliptical, NGC 3550, 
to be 5ai - 6.5, which is in the normal range for a galaxy of 
its luminosity. For the other galaxies, none of whose globular 
cluster systems have been studied to date, we used the empir- 
ical relation between specific frequency and V-band galaxy 
luminosity determined by Peng et al. l2008l from observations 
of 100 early-type galaxies in the ACS Virgo Cluster Survey 
(Cote et. al'2003^, with S n values ranging from ^1-5 and 
the V-band magnitudes measured directly from the CFH12k 
images of Andreon et al. 120071 using SExtractor This rela- 
tion agrees well with that found for other galaxies (see, e.g.. 
Ashman & Zepf l 19981 and references therein) and hence pro- 
vides a reliable input for the simulations. A total number 
of globular clusters was assigned to each galaxy in Table 1 
based on its luminosity and assumed specific frequency. 

Each simulated globular cluster was assigned a luminos- 
ity by randomly sampling a Gaussian luminosity function 
with a mean absolute magnitude of Mj - -8.4 mag and a 
standard deviation cr = 1.4 mag, which is also consistent 
with previous observations of the GC population of NGC 
3550 (Blakeslee et al. 119971 1. Clusters whose apparent mag- 
nitude at the distance of A1185 would be fainter than our 
Ifsi4w = 27.3 mag completeness limit were then removed 
from further consideration. Although Jordan et al. 120071 and 
Villegas et al. 2010 have shown that the properties of the 
GCLF correlate with parent galaxy luminosity, it is nearly 
uniform among high-mass galaxies that contribute the bulk 



of contaminating globular clusters, with turnover magnitudes 
and dispersions in agreement with the values adopted here. 
The five most luminous galaxies in Table 1, for example, 
which contribute the lion's share of halo globular clusters, 
would be expected to have cr ranging from ~ 1.2 to 1.4 mag. 
Hence, any non-universality of the GCLF can safely be ig- 
nored in the present analysis. 

The spatial distribution of globular clusters in galaxy halos 
varies widely but is generally more extended than that of the 
starUgh t (e.g Harris [19911 Ashman & Zepf 1998 ; Brodie & 
S trader I2OO6I 1. We assumed that the globular cluster system 
around each galaxy has a spherical distribution that follows 
a Sersic d 19681 1 profile 



I(r) = I, exp[-b„(r/re) 



I/n . 



1] 



(3) 



where is the projected surface brightness at the effective 
radius, , that contains half the total light of the system, n is 
the Sersic index that determines the overall shape of the pro- 
file, and b„ is a parameter related to n (see Graham & Driver 
,2005 and references therein). Surface brightness profiles for 
the stellar distribution in each galaxy were measured on the 
CFH12k R-band image of A1185 (Fig.[T]l using the ellipse 
routine in IRAF, and Sersic parameters were then obtained 
by fitting the profiles to equation 3 using minimization}^ 
As an initial guess, the simulations assumed that the effective 
radius of each galaxy's globular cluster distribution is twice 
that of its stellar distribution (Peng et al. 2011, in prepara- 
tion). 

Tidal forces from the gravitational field of A1185 and 
galaxy-galaxy interactions are likely to limit the radial extent 
of each galaxy's globular cluster distribution (e.g., Merritt 
1984 although see Mihos et al. 120051 and Janowiecki et al. 
12010) . For this reason we simulated three different scenar- 
ios. For one, the globular cluster distribution was not trun- 
cated, which likely leads to an overestimate of the amount of 
contamination from halo globular clusters in our ACS field. 
For the other cases, the globular cluster population of each 
galaxy was truncated beyond a radius of 50 kpc or 100 kpc. 
Previous observations by Blakeslee d 199 61) showed that the 
number density of point sources surrounding NGC 3550 be- 
comes constant beyond a projected distance of 100", sug- 
gesting that this galaxy's globular cluster population does not 
extend beyond ~ 63 kpc. Similarly, in a wide-field imaging 
survey Rhode & Zepf (2001 ; 2004 ) found that the maximum 
radial extent of the globular cluster systems of four early- 
type galaxies in the Virgo cluster ranged from ~ 30-100 
kpc. The assumed radial distribution of globular clusters in 
galaxy halos is the greatest uncertainty in these models. 

Using the procedures and input parameters described above, 100 
random realizations of the expected distribution of globular clus- 
ters in the halos of the galaxies in Table 1 were generated and 



the simulated clusters brighter than If^i 



27.3 mag that fall 



within the ACS field boundaries were counted. An example is 
shown in Fig. 1 1 

Based on these simulations, we estimate that ~ 711+26 
IGC candidates in the Al 185 core field are likely to be globular 



As Fig.[lO|shows, NGC 3550 is actually a merging system of three 
galaxies within a common envelope, a brightest cluster galaxy caught 
in the act of formation. The light from the two companion galaxies was 
subtracted before measuring the underlying surface brightness profile 
of NGC 3550. 
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Table 1. Galaxies near the center of A 11 85 used to estimate contamination in the ACS field from halo globular clusters 



Galaxy 


RA (2000) 


DEC (2000) 


V (mag) 


z 


NGC 3550 


11:10:38.7 


28:46:06 


13.2 


0.0346 


NGC 3552 


11:10:42.8 


28:41:35 


14.3 


0.0331 


CGCG 155-081 


11:10:37.5 


28:43:58 


14.9 


0.0310 


MCG+05-27-004 


11:10:52.0 


28:42:09 


15.3 


0.0340 


NGC 3553 


11:10:40.5 


28:41:05 


15.3 


0.0326 


2MASX Jl 1 105941+2844424 


11:10:59.4 


28:44:43 


15.7 


0.0327 


2MASX Jl 1 104649+2841338 


11:10:46.5 


28:41:34 


15.9 


0.0326 


2MASX Jl 1 103987+2842248 


11:10:39.9 


28:42:25 


16.3 


0.0288 


2MASX Jl 1 103964+2845208 


11:10:39.6 


28:45:21 


16.9 


0.0383 


2MASX Jl 1 103683+2842138 


11:10:36.8 


28:42:14 


16.9 




2MASX Jl 1 104185+2845408 


11:10:41.9 


28:45:41 


17.0 




2MASX Jl 1 105850+2842384 


11:10:58.5 


28:42:39 


17.1 


0.0346 


SDSS Jl 1 1035.04+284205.2 


11:10:35.0 


28:42:05 


17.3 


0.0338 


2MASXJ1 1103767+2841528 


11:10:37.7 


28:41:53 


17.6 


0.0343 



clusters in galaxy halos if their spatial extent in unlimited, most 
of them belonging to NGC 3550. We emphasize, however, that 
this probably overestimates the number of contaminating halo 
globular clusters. For the case in which each galaxy's globular 
cluster system is assumed to be tidally truncated beyond 100 
kpc, the estimated number of contaminating objects is ~ 472 ± 
30, and for a limiting radius of 50 kpc the number falls to ~ 
166 ± 14. 

Combining the estimated contamination from globular clus- 
ters in galaxy halos and the expected number of unassociated 
objects along the line of sight, we estimate that less than half of 
the 1422 IGC candidates brighter than If&uw - 27.3 mag can be 
accounted for this way, and hence the majority of detected ob- 
jects are likely to be genuine IGCs. Our best estimate, based on 
the assumption that globular clusters are not found beyond 100 
kpc in galaxy halos, is that ~ 1616 ± 91 IGCs of all magnitudes 
exist in the core ACS fieldE] 

However, although there is no doubt that we have detected a 
large population of globular clusters in the core of A1185, it is 
diflicult to say with certainty whether any of them are truly inter- 
galactic. In addition to the aforementioned assumptions and un- 
certainties inherent in our simulations, other factors might also 
plausibly raise the number of globular clusters associated with 
galaxies. For example, although low- or intermediate-luminosity 
galaxies (fainter than My - -18) have not been included in the 
simulations because their globular cluster populations are ex- 
pected to be sparse, a large population of such galaxies could 
collectively contribute several hundred additional contaminating 
halo clusters scattered over the ACS field. Uncertainties in the 
magnitudes of the galaxies in Table 1 could also raise or lower 
the expected number of galactic globular clusters. 

It is, of course, possible to boost the predicted number of 
contaminating objects to match the number of detected objects 
in the ACS field by arbitrarily increasing the globular cluster 
populations of the galaxies or changing their radial distributions. 
For example, doubling the assumed specific frequency of every 



galaxy in Table 1 and assuming that assuming that the globu- 
lar cluster distribution is twice as extended as the galaxy light 
yields an estimated 1292 + 34 contaminating halo globular clus- 
ters in the ACS field. However, such a model seems unlikely, 
especially in light of the measured specific frequency of NGC 
3550 (Blakeslee et al. 1997). To further assess the validity of 
such a high specific frequency model, the simulated distributions 
of globular clusters (including the addition of 142 ± 12 randomly 
distributed foreground/background objects) were compared to 
the observed distribution using a two-dimensional Kolmogorov- 
Smirnov test (Peacock 119831 Fasano & Franceschini 119871 ) as 
implemented by Press et al (1992). Based on this analysis, the 
probability that the simulated and observed globular cluster dis- 
tributions could both be drawn from the same parent distribution 
is much less than 1%, although one must bear in mind the limi- 
tations of these models. 



Given the uncertainties in the simulations, and to explore the 
full range of parameter space - including more extreme possibil- 
ities - several thousand additional simulations were run in which 
the parameters of each galaxy's globular cluster system were 
randomly chosen in the range 1 < S n < 20, 0" < < 200" 
and < n < 10. These experiments demonstrate that, even in the 
most favorable case in which the spatial extent of each galaxy's 
GC system is not tidally limited, a minimum necessary condition 
for ~ 1000 - 1500 contaminating halo globular clusters to fall 
within our central ACS field is that two or more of the bright- 
est member galaxies (NGC 3550, NGC 3552, CGCG 155-081 
and MCG +05-27-004) would have to possess extremely rich 
globular cluster systems, with S n ^ 10-20 (see Fig. 12 for an 
illustration). While this possibility cannot be ruled out, it seems 
unlikely that Al 185 would have the distinction of being the only 
known cluster with multiple high S n galaxies, especially given 
the unexceptional size of NGC 3550's globular cluster popula- 
tion (Blakeslee et al. |1997l ). 



^ This estimate of ~ 1616 IGCs is based on the detection of 1422 
objects down to Ifhaw - 27.3 mag, subtracting 142 background ob- 
jects brighter than this same limiting magnitude and assuming that 472 
of the IGC candidates are globular clusters in galaxy halos. The re- 
sulting number is then doubled to account for the fact that we observe 
only the brighter half of the GCLF. Formal statistical uncertainties in 
the number of detected IGCs, unrelated point sources, and halo glob- 
ular clusters have been added in quadrature, however systematic errors 
resulting from our lack of knowledge about the globular cluster systems 
of individual galaxies make the results considerably more uncertain. 



We conclude that, for any reasonable assumptions, it is dif- 
ficult to account for the large number of globular clusters de- 
tected in our central ACS field without including a significant 
population of IGCs that do not belong to galaxies. Although we 
cannot completely iTile out the possibility that most or all of the 
globular clusters detected in the core of Al 185 could be associ- 
ated with galaxies, this seems unlikely. Further insights into the 
nature of these objects could come from their luminosities and 
colors, which are examined in more detail below. 
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Fig. 11. Simulations of the expected contamination from globular clusters in halos of galaxies that lie near the A1185 core field. The 
left panels shows the real distribution of globular clusters brighter than If&i4w - 27.3 mag in the ACS image of the heart of 1 185. 
The right panel shows a superposition of 10 simulated distributions of halo globular clusters associated with the galaxies listed in 
Table 1 down to the same If&hw - 27.3 mag limit and assuming no truncation of their radial distribution, as described in the text. 



3.2. The globular cluster luminosity function 

The luminosity function (or, equivalently, the mass function) is 
a fundamental property of globular cluster systems. The present- 
day GCLF has been shaped by many factors, including physical 
conditions at the time of their formation, internal processes such 
as stellar mass loss and two-body relaxation, and stripping or 
destruction of clusters in galaxy environments as a result of disk 
shocking, tidal forces, and dynamical friction. Given this variety 
of influences, it is remarkable that the GCLF in galaxies can, to 
first order, be described universally by a Gaussian distribution 
with a maximum at Mj ^ -8.4 and a dispersion of cr ^ 1.2-1.4 
magnitudes (e.g., Harris fTWTT l. However, some correlation is ob- 
served between GCLF properties and parent galaxy luminosity 
(Jordan et al. l2006l 120071 ViUegas et al.'SOTO'), which suggests 
that environment plays a role in the formation and/or evolution 
of globular cluster systems. Several authors have shown, for ex- 
ample, that an initial power-law mass function can be trans- 
formed into one with a log-normal shape as a consequence of 
cluster disruption and evaporation in galaxies (see, e.g.. Fall & 
Zhang |200 1; Ves perini et al. l2003l McLaughlin & Fall l2008l 
Elmegreen EoTOb . 

It therefore seems plausible that the luminosity function of 
IGCs born outside of galaxies might diff'er from that of globular 
clusters which formed and evolved within galaxies, although it 
is hard to predict a priori what form such differences might take. 
If, on the other hand, IGCs originated in galaxies and were sub- 
sequently stripped from halos or spilled into intergalactic space 
during the disruption of their parent galaxies then the IGC lumi- 
nosity function is likely to be similar to that of galactic globular 
clusters. Hence the luminosity function of globular clusters in 
the core of A 11 85 might provide some insights into their origin. 



Figure 13 shows the observed luminosity function of glob- 
ular cluster candidates based on the deep F814W observations. 
The number of objects in each luminosity bin has been corrected 
by subtracting the expected contamination from foreground stars 
and unresolved backg^und galaxies (Fig. [8]l and correcting for 
incompleteness (Fig. M. Using minimization to fit the ob- 
served GCLF to a Gaussian distribution with the turnover mag- 
nitude fixed at IfsiAw = 27.3 mag but with the amplitude and 
width of the distribution as free parameters yields a best fit with 
cr ^ 1.2 ± 0.2, which is shown for comparison in the figure. 
Within the uncertainties, this falls well within the normal range 
for GCLFs in galaxies. 



3.3. The globular cluster color distribution 

It is well established that the colors of old globular clusters, 
which are determined primarily by their metallicities, show a 
dependence on parent galaxy luminosity, in the sense that lumi- 
nous galaxies have redder (more metal-rich) globular clusters on 
average than fainter galaxies (e.g., van den Bergh 1975; Brodie 
& Huchra [TWll Forbes & Forte 120011 Cote, Marzke & West 
|I998; Peng et al. 120061 1. Hence the color distribution of glob- 
ular clusters in the core of A1185 could, in principle, provide 
some constraints on the numbers and types of donor galaxies 
from which an IGC population might have been stripped (e.g., 
Forbes, Brodie & Grillmair 1997; Cote, M arzke & West 1991 
Larse n et al. 2001; Gebhardt & Kissler-Patig [T999l Coenda et al. 
120091 ). In reality, however, such constraints are likely to be weak 
because of degenerate dependencies of globular cluster colors on 
galactocentric distances and host galaxy properties. A prepon- 
derance of blue IGCs could arise, for example, if the majority 
came from dwarf galaxies that were tidally disrupted or if they 
were stripped from the outer halos of giant galaxies whose glob- 
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Fig. 10. NGC 3550, the brightest member galaxy in Al 185. The 
image on top is displayed with a linear intensity scale and NGC 
3550 appears as a single large elliptical galaxy. The logarithmic 
display below reveals that NGC 3550 is actually a complex sys- 
tem of three smaller galaxies within a common envelope. 



ular clusters are bluer on average than their counterparts closer 
to the galaxy center (Geisler et al. [19961 Rhode & Zepf l2004l 
Bassino et al. 120061 Tamura et al. I2006I I. Nevertheless, if the 
globular clusters are found to have a distinctive color distribu- 
tion then this could be an important clue to their origin. 

Artificial star experiments show that matched detection of 
objects added to both the F555W and F814W images is 100% 
complete to Ifs\4w - 26.5 mag, as the shallower F555W ob- 
servations result in fewer detections at fainter magnitudes. To 
minimize the effect of this incompleteness on the inferred color 
distribution of globular clusters, we consider only those objects 
brighter than this magnitude limit. 

Figure [14] shows the observed color distribution of globular 
clusters in the Al 185 core field. These colors are quite similar to 
those of globular clusters associated with galaxies (Gebhardt & 
Kissler-PatiglHSSl Peng et al l2006l Forbes et al. l2U07] l. skewed 
strongly toward bluer (more metal-poor) colors. The peak in the 
color distribution occurs at V - / ^ 1.05, in excellent agree- 
ment with the results of Gebhardt & Kissler-Patig 119991 who 
found peak colors around V - 1 - 0.94 -1.13 depending on host 
galaxy luminosity. Although it is impossible to distinguish be- 
tween IGCs and halo globular clusters in the HST images, it is 
clear that if IGCs are present then their colors cannot be grossly 
different from those of normal globular clusters. The colors are 
predominantly blue, as expected if IGCs are the surviving or- 



relative Sfj 



Fig. 12. An illustration of the expected contamination from glob- 
ular clusters in galaxy halos near the center of Al 185 and its de- 
pendence on the assumed properties of those systems. The glob- 
ular cluster population of each galaxy in Table 1 was simulated 
by assuming that its radial distribution is described by the same 
Sersic index n as the galaxy's stellar light, but with an effective 
radius that is some multiple of the galaxy's (the axis labelled 
'relative r/). Similarly, the number of globular clusters belong- 
ing to each galaxy was allowed to vary relative to the expected 
number if the galaxy had a normal specific frequency for its lu- 
minosity as summarized in Peng et al. 120081 (the axis labelled 
'relative S n')- Hence, for example, a relative - 2 means that 
the effective radius of a galaxy's globular cluster distribution was 
assumed to be twice that of its stellar distribution, and an 'effec- 
tive S N = 3' means that a galaxy was assigned three times the 
'normal' number of globular clusters for a galaxy of its lumi- 
nosity. Additional details of the simulations are discussed in the 
text. The third axis shows the average number of contaminating 
globular clusters brighter than Ifs\4w ^ 27.3 mag predicted to 
fall within the central ACS field in Al 185 based on many simu- 
lations. Although the relative and S n were allowed to vary in- 
dependently for each galaxy, for clarity the results shown in this 
figure assume constant values of relative and relative S n for 
all galaxies. These simulations demonstrate that the large num- 
ber of globular clusters seen in our ACS field cannot all reside in 
galaxy halos unless most or all of the galaxies in Table 1 have 
unusually rich globular cluster systems. 



phans of disrupted dwarf galaxies or were stripped from the 
outer halos of larger galaxies. 



Figure 14 hints at possible bimodality or asymmetry of the 



globular cluster color distribution. To quantify this impression, 
we used the KMM mixture mode algorithm (McLachlan & 
Basford '1988) as implemented by Ashman, Bird & Zepf ll994l 
KMM uses a maximum likelihood estimator to determine the 
best fitting mixture of one or more Gaussian components to an 
observed distribution and estimates the mean and fraction of data 
points belonging to each component. Because outliers can ad- 
versely affect the KMM results (see Ashman et al. 1994'for a dis- 
cussion) we restricted our analysis to the 616 clusters with colors 
in the range 0.8 < V-I < 1.4 and /f814H' ^ 26.5 mag. Assuming 
homoscedasticity (i.e., the same variance for each subpopula- 
tion), the KMM algorithm strongly rejects (at a confidence level 
of more than 99%) the hypothesis that the color distribution seen 
in Fig. 14 is adequately described by a single Gaussian. The best 
fitting model has two Gaussian components, with roughly 2/3 of 
globular clusters belonging to a blue subpopulation with a mean 
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Fig. 13. The observed luminosity function of globular clusters 
in the core of A 11 85. The data have been corrected for in- 
completeness and the expected number of contaminating fore- 
ground/background objects have been removed based on the av- 
eraged counts at each magnitude from the four control fields. 
The luminosity function data are shown only to If&uw - 27.8 
mag because beyond that incompletenss corrections become too 
large (more than a factor of two) to be reliable. The red curve is 
the best-fitting Gaussian luminosity function with a peak lumi- 
nosity at /f8i4H' - 27.3 mag and a standard deviation cr - 1.2, as 
discussed in the text. The observed luminosity function is clearly 
consistent with the standard GCLF observed in many galaxies. 
The inset shows the raw luminosity function of globular cluster 
candidates without any corrections for background contamina- 
tion or incompleteness. 



color (V-I) - 1 .04 and 1/3 belonging to a redder subpopulation 
with<y-/) ^ 1.23. 

Given the detection of distinct blue and red subpopulations, 
we looked for evidence of the so-called 'blue tilt' in the color- 
magnitude relation (Fig.|9]l, an observed tendency for blue glob- 
ular clusters to become progressively redder at brighter magni- 
tudes (e.g., Haiiis et al. 120061 Stra der et al. |2006; Pen g et al . 
12009 '; Blakeslee, Cantiello & Peng l2010l Mieske et al. l20T0l l. 
In principle, the strength of the blue tilt, which is usually inter- 
preted as a mass-metallicity relation, could provide clues to the 
origin of the globular clusters in the center of Al 185. As noted 
by Mieske et al. (2010), for example, the slope of the blue tilt in 
the Virgo and Fornax clusters of galaxies varies as a function of 
host galaxy luminosity and galactocentric distance. 

For this purpose, we divided the sample of 616 globular clus- 
ters satisfying the aforementioned color and magnitude criteria 
into four luminosity bins containing equal numbers of objects. 
The KMM algorithm was then used to find the best fitting blue 
and red peaks in the color distribution for each luminosity bin 
(again assuming homoscedasticity) with V - I - 1 .04 and 1 .23 
chosen as the initial guess for their locations. Within each lumi- 
nosity bin, the color distribution is well fitted by a mixture of two 
Gaussians. Because the KMM analysis does not provide a direct 
estimate of the uncertainties in the derived quantities, bootstrap 



Fig. 14. The observed V-I color distribution of globular clusters 
in the core of A1185. The bar indicates typical colors of glob- 
ular clusters in the sample of Gebhardt & Kissler-Patig ( 1999l l, 
corresponding to the FWHM of their color distribution). 



resampling was used to estimate the variance in the mean colors 
of the blue and red subpopulations. Results are shown in Fig. 15 
Despite the relatively weak sensitivity of V - / colors to metal- 
licity (compared to the longer wavelength baselines provided by 
g - ZOT B - I used in most previous studies) a blue tilt is appar- 
ent. No obvious trend is seen for the red clusters. These visual 
impressions are confirmed by statistical analysis, which yields a 
Pearson correlation coefficient r - -0.99 (highly significant) for 
the color-magnitude relation of the blue subpopulation but only 
r = -0.45 (not significant) for the red subpopulation. The slope 



of the best fitting linear regression to the blue tilt seen in Fig. 15 
is d{V - I)/dI ^ -0.039 ± 0.013, with the uncertainty estimated 
from bootstrap resampling. This slope is within the range found 
in previous studies. Given the sizable uncertainties, however, we 
refrain from further analysis or interpretation of these results. 

4. Discussion 

The results presented here show that the ~ 1300 globular clus- 
ters found in our ACS field are indistinguishable, at least in 
terms of their optical colors and luminosities, from the metal- 
poor globular cluster populations observed in hundreds of galax- 
ies to date. Although a non-negligible fraction of these objects 
are undoubtedly chance projections of globular clusters resid- 
ing in the extended halos of galaxies, we have argued based on 
detailed simulations that the majority are likely to be intergalac- 
tic. Nevertheless, given the uncertainties in these simulations we 
cannot confidently rule out the possibility that most or even all 
of the detected globular clusters might belong to galaxies rather 
than roaming the space between them. 

Confirmation that the candidate IGCs are truly intergalactic 
could be obtained in several ways. The most direct, but observa- 
tionally challenging, method would be to measure velocities of 
these objects to determine whether they are indeed gravitation- 
ally bound to individual galaxies or whether they constitute an 
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Fig. 15. Color-magnitude diagram for 616 globular clusters in 
the core of A1185 within the range 0.8 < V - / < 1.4 and 
/f814w ^ 26.5 mag. The blue and red points indicate the lo- 
cations of the blue and red peaks determined by KMM fitting 
of the color distribution in four different luminosity bins to a 
model with two Gaussian components, as described in the text. 
Uncertainties in the peak locations were calculated via bootstrap 
resampling. The dashed lines denote the locations of the blue 
and red peaks found from KMM analysis for the entire sample 
of 616 clusters. A correlation between color and magnitude is 
apparent for the blue globular clusters (dubbed the 'blue tilt') 
while no such trend is seen for the red clusters. 



independent dynamical system in the core of All 85. In particu- 
lar, the largest potential source of contaminating halo globulars, 
NGC 3550, has an observed radial velocity of 10,520 km s"' 
(Tonry [1985] gives velocities of 10447, 10388 and 1100 km s"' 
for each of the three components seen in Fig. 9), which differs 
by ~ 1200 km s"' from the 9300 km s"' mean velocity of the 
A1185 system (Mahdavi et al. 1996), consistent with the spatial 
and kinematic substructure seen in this dynamically young clus- 
ter Consequently, halo globular clusters associated with NGC 
3550 could be readily distinguished from genuine IGCs moving 
in the gravitational potential of A1185 on the basis of their ex- 
pected large velocity differences. The faint magnitudes of even 
the brightest IGC candidates in All 85 make such observations 
challenging, though not impossible, with 8-m to 10-m class tele- 
scopes. 

Deep images of additional fields surrounding NGC 3550 
could also help to determine whether IGCs are present. Globular 
clusters that are gravitationally bound to galaxies should gener- 
ally be distributed symmetrically about their centers. If the num- 
ber density of globular clusters in the central Al 185 field shown 
in Fig. [T| is found to be significantly higher than that in other 
fields equidistant from NGC 3550, then this would support the 
idea that they are intergalactic. Several strategically placed fields 
might even allow the centroid of the putative IGC population to 
be localized. 

Additional clues to the nature and origin of the globular clus- 
ter population in the heart of A1185 could come from mea- 
surements of their ages and metallicities. Although it is difficult 



at present to obtain such information spectroscopically because 
of the distance to A1185, a combination of optical and near- 
infrared colors can provide an alternative measure of ages and 
metallicities (e..g, Puzia et al. 120021 Hempel et al. 120071 Carter 
etal. 120091) . 

Finally, we note that further study of Al 185 may offer a tan- 
talizing opportunity to glimpse the ongoing creation of IGCs. 
The merging galaxy system Arp 105 (see Fig. [T} is ejecting 
a long, narrow plume of material that probably includes ex- 
tant globular clusters destined to become future intergalactic 
vagabonds, as well as newly created dwarf galaxies and pos- 
sibly young star clusters (Due & Mirabel 1994). Similar tidal 
streams in other nearby galaxy clusters (e.g., Gregg & West 
[T9981 Calcaneo-Roldan et al. l2550l Rudick et al. l2009T l would be 
excellent places to search for other IGCs. Deep, large-area sur- 
veys such as the ongoing Next Generation Virgo Cluster Survey 
(Ferrarese et al. 2010, in preparation) and the ACS Treasury 
Survey of the Coma Cluster (Carter et al. 2008) will also pro- 
vide rich databases that can be searched for IGCs. 



5. Conclusions 

Using deep HST/ACS images, we have confirmed the presence of 
several thousand globular clusters in the heart of the rich galaxy 
cluster Al 185, in a field that is not centered on any bright galaxy. 
Simulations suggest that the majority of these objects may be 
intergalactic, although the large uncertainties in these models 
make it impossible to reach firm conclusions about the reality 
or size of such a population. Our best estimate is that a total of 
~ 1600 IGCs reside in the core ACS field, although we cannot 
rule out the possibility that they could be gravitationally bound 
to galaxies. The colors and luminosities of these objects are sim- 
ilar to those of metal-poor globular clusters found in galaxies, 
suggesting that - if they are indeed intergalactic - they proba- 
bly formed in galaxies and were removed later by tidal stripping 
or other dynamical processes. We suggest that future observa- 
tions, especially velocity measurements, could establish the true 
nature of the globular clusters in the center of A1185 as inter- 
galactic wanders or galactic residents. 
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